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<N : Abstract 

^ , Populations of swimming microorganisms produce fluid motions that lead to dramatically enhanced dif- 
fusion of tracer particles. Using simulations of suspensions of swimming particles in a periodic domain, we 

■ capture this effect and show that it depends qualitatively on the mode of swimming: swimmers "pushed" 
q— . from behind by their flagella show greater enhancement than swimmers that are "pulled" from the front. 

§ I The difference is manifested by an increase, that only occurs for pushers, of the diffusivity of passive tracers 

■ and the velocity correlation length with the size of the periodic domain. A physical argument supported by 

B' 

I . a mean field theory sheds light on the origin of these effects. 
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In the present work we consider the possibility of emergence of large scale fluid motion and 
enhanced fluid transport in populations of small swimming organisms. At the global scale, it has 
been suggested that swimming organisms such as krill can alter mixing in the oceans [fl, Q]. At 
the laboratory scale, experiments with suspensions of swimming cells have revealed characteristic 
swirls andiets much larger than a single cell, as well as causing increased diffusivity of tracer 
particles [3, 4, si [5, Q]- This enhanced diffusivity may have important consequences for how 
cells reach nutrients, as it indicates that the very act of swimming towards nutrients alters their 
distribution. The enhanced diffusivity has also been proposed as a scheme to improve transport in 
microfluidic devices [[st] and might be exploited in microfluidic cell culture of motile organisms or 
cells. 

The feedback between the motion of swimming particles and the fluid flow generated by that 
motion is thus very important. Nevertheless, in the literature on collective dynamics of self- 
propelled particles, this effect has received little attention. Most previous attempts to understand 
the collective motion of swimming micro-organisms fall into two broad categories: "particle level" 
approaches based on ad hoc interaction rules between the individual agents, e^.J9(]j_and contin- 



uum models of active suspensions based on phenomenological field equations IIO , 



11 



121] . The hy- 



drodynamics of single swimming organisms at low Reynolds number have been studied for many 
years 113, 14, isl, and studies of hydrodynamic interactions (HI) have also been performed for 



pairs of swimmers 



OIlV, 



1811 . However, only a small number of studies ll 
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have examined particle level models of populations of swimmers that include the fluid motion 
caused by the swimmers and its influence on transport. Hernandez-Ortiz et al. Iioi I20I1 devel- 
oped a simple physical model of self-propelled particles, and performed simulations of them in 
a confined domain. Their work showed that multi-body HI between self-propelled particles were 
sufficient to lead to fluid motions characteristic of those seen in experiments. Llopis and Pag- 
onabarraga [12 ill performed simulations of self-propelled particles in two-dimensions, observing 
aggregates of swimming particles and a decrease of velocity with increasing volume fraction. 
Saintillan and Shelley 12211 performed simulations of long, slender self-propelled rods propelled 
by an imposed shear stress on the surface of part of the rods. Their main result was that initially 
nematically ordered swimmers do not remain aligned except with very near neighbors. Pedley and 
coworkers 11231 12411 examined the rheology and random walk motion of semi-dilute suspensions of 
self-propelled spheres. 

We consider here a suspension of neutrally buoyant rodlike swimmers in a spatially periodic 
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FIG. 1 : (color online) Illustration of a pushing organism, swimmer model, and fluid disturbance they cause. 
The double-arrow signifies the flagellum force acting on the bead and opposite force acting on the fluid, 
both acting at the center of the first bead. The blue lines represent streamlines of the axisymmetric fluid 
disturbance. A puller produces the same streamlines with the arrows reversed. 



cubic fluid domain of side length L, in a concentration regime dilute enough that the dominant 
interactions between the particles arise from the fluid motions that they generate as they swim. 
Each swimmer has a characteristic length i, and in isolation moves in a straight line with speed 
Vis- It is assumed that the Reynolds number, VisL/u, is much less than 1, where u is the fluid 
kinematic viscosity. The concentration is given as an effective volume fraction 0e = uNi^/ {6L^) 
- this would be the true volume fraction if the swimmers were spheres of diameter i. To allow 
large populations (> 10'' swimmers) over long times, a simple model of each swimmer is adopted, 
following our previous work lll9ll . 

Each self-propelled particle is represented as two beads connected by a stiff spring with equi- 
librium length i as shown in Figure 1 . The unit vector pointing along the axis from bead 1 to 
bead 2 is denoted n. The propulsion is provided by a "phantom flagellum" that we do not treat 
explicitly, but only through its effect on the swimmer body and the fluid. This flagellum exerts a 
force F/ on bead 1 of the swimmer, and also exerts a force —Fj on the fluid. This force on the 
fluid occurs at the position of bead 1 . With this model we can consider "pushers" or "pullers" 
depending on whether Fj is parallel or anti-parallel to n, respectively. A pusher sends fluid away 
from it fore and aft, with fluid moving toward its "waist", and vice versa for a puller. A cell whose 



flagella propel it forward predominantly from behind would be a pusher il9 



251]. We shaU see that 



this distinction is important. The swimmers also interact with an excluded volume potential using 
the repulsive portion of the Gay-Berne potential 1261] . 
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The motion of the swimmers is calculated by writing a force balance (neglecting inertia because 
of the small size of a microorganism) for each bead. On bead 1 of each swimmer this is F/ + Fhi + 
Fci +Fei = 0, where F^i is the hydrodynamic drag force on the bead, F^ is the connector (spring) 
force on the bead, and Fei is the excluded volume force on the bead. On bead 2 of each swimmer 
the force balance is identical except without a flagellum force. The drag force is written using 
Stokes law, Fhi = —6nT]a[ii — v'(ri)], where a is the bead radius, r] is the fluid viscosity , ri 
is the position of the bead, and v'(ri) is the fluid velocity at ri generated by all other beads and 
phantom flagella. This fluid velocity is calculated using an order N method [|27ll . treating each 
bead of the swimmer as a regularized point force. We set i = 3a. Note that the net force exerted 
by an isolated swimmer on the fluid is zero - overall, a neutrally buoyant swimmer is a force 
dipole to leading order. This model captures the universal far-field behavior while neglecting the 
near-field corrections to interactions between swimmers, which are dependent on the details of the 
organism. The validity of this approximation is supported by recent simulations ll28n . Finally, we 
note that the limited set of results with multi-bead rod swimmers is qualitatively consistent with 
those for the two-bead swimmers. 

The fluid motion generated by each swimmer perturbs the trajectories of fluid elements (tracers) 
and other swimmers. At the volume fractions considered here, this motion also eliminates any long 
range orientational order among the swimmers. The motion of both tracers and swimmers, while 
ballistic at short times, becomes diffusive at long times. The tracer diffusion in particular shows 
important and unanticipated behavior. Figure 2 shows tracer mean squared displacements (inset) 
and the long-time tracer diffusivity Dt in a suspension of pushers as a function of system size 
for (pe = 0.1. The diffusivity grows significantly with system size, greatly enhancing it beyond 
that for a suspension of pullers (also shown), which exhibits only weak or negligible system size 
dependence. 

The system size dependence observed for pushers is due to collective behavior; simulations 
in which the swimmers do not interact, swimming in straight paths, do not show a signifi- 
cantly changing diffusivity. This observation has important implications for transport of nutrients, 
chemoattractants, and other chemical species in the environment of a population of swimming 
organisms. For example, organisms which concentrate toward an attractant will disperse that at- 
tractant more rapidly if they are pushers than if they are pullers. The remainder of this report 
describes efforts to gain an understanding of the increased transport. 



The diffusivity is related to the fluid velocity through the Green-Kubo relation [|291],Dt 
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FIG. 2: (color online) System size dependence at volume fraction <j)e = 0.1 for pushers. Variation of 
mean squared tracer velocity {v^) (squares), [Vg) — vf^ (triangles), tracer diffusivity Dt (circles), and tracer 
correlation time (stars). The dashed line is a power law with a = 0.63. The tracer diffusivity in a suspension 
of pullers is shown with diamonds, (inset) Tracer mean-squared displacements versus time for independent 
swimmers (blue), pullers (orange), and pushers (green) for N = 400 (solid) and = 3200 (dashed). 
Averages over 400 tracers typically give statistically converged results. 

(1/3) /Q°°(vf(0) ■ \i{t))dt, where is the velocity of a tracer (fluid element) and angle brackets 
indicate an ensemble average. The mean squared tracer velocity {vf) = (vt(0) ■ Vf(0)) is shown in 
Fig. 2 and has a much weaker system size dependence. The simple idea that the swimmer velocity 
is the isolated value plus the tracer velocity, = Vj^ + , predicts that (f^) ^ vf^ + {v^), which 
is verified in Fig. 2. 

We now define a correlation time, r = Dt/ (v^), for the trajectories of tracer particles - tracer 
motion becomes diffusive on times t ^ r. Figure 2 shows r versus L; over the entire range 
of system sizes considered, this relation displays a clear power law dependence r ~ L"", with 
a = 0.63 ± 0.03, over a range spanning a factor of 8 in L and thus 512 in A^. This data represents 
the first examination of system size effects on the fluid response in swimming suspensions, and 
includes orders of magnitude more swimmers than any previous studies. For larger systems the 
behavior may transition to a non power law regime, but the system size dependence shown here 
will nevertheless cause significant enhancement of transport in a suspension of pushers. 

The time r represents the time a tracer moves in the same direction. A corresponding step 
length or correlation length It for the tracer's random walk scales as T{vfy/'^. Because (vf) has a 
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much weaker dependence on L than r, our resuks imply the existence of a correlation length for 
pushers that obeys U ~ t{L/ tj"' ~ ^(i-")/,^^ over the range of system sizes examined. The suspen- 
sion of pushers leads to enhanced correlation lengths in the fluid, beyond pullers or independent 
swimmers, because of this increase with system size. 

Before elaborating on correlation lengths, we briefly address the concentration dependence, 
using simple scaling arguments. In the Green-Kubo relation, the tracer velocity is a sum of dis- 
turbances due to each swimmer. Assuming independence of the swimmers in the dilute limit, the 
cross terms from the dot product vanish, resulting in Dt ~ 0e. For the same reason, {v^) ~ 0e> 
which leads to a correlation time r that is independent of volume fraction. Turning to swimmers 
rather than tracers, we note that in the dilute limit, the swimmers continue in straight lines until 
perturbed by other swimmers. If we define a cross section a for these "collisions", the mean free 
path or step length Ig scales as £^/(cr0e)- Because l^vf) vf^ + {vj) and i^v^) ~ 0e, the root- 
mean-squared swimmer velocity is (f^)^^^ ~ Vis + 0{(j)e)- Therefore, the diffusivity scales as 
Da ~ {vlY^'^ls ~ For 0e ^ 10^^ simulation results follow these predictions (not shown), 
while above that volume fraction, they begin to deviate, indicating a breakdown of the assumption 
of independent swimmers. 

We now return to the issue of correlation lengths, and consider the fluid velocity spatial cor- 
relation function, C(r) = {L^"^ J vt(r') • vt(r' + r)(ir'). This is shown in Figure 3 at a volume 
fraction of (p,. = 0.1 for a range of system sizes. For independently distributed point dipoles, this 
correlation is Cjd(r) = Vji/r, where = (peCp / (Sr^^Tr^f"^), and d is the dipole moment 1301] . Be- 
cause of the finite size of the swimmers and domain, real independent swimmers will deviate from 
this relation when r < £ or r ^ L. Even this simple result illustrates the possible complexities that 
arise in considering swimming suspensions; the velocity correlations are very long -ranged even in 
the absence of correlations between the positions and orientations of the swimmers. 

Figure 3 shows the correlation function at a characteristic system size of = 3200 (L = 
25.589£), for pushers, pullers, noninteracting pushers (where the swimmers do not feel the fluid 
motion but the tracer particles do) and the point dipole theory result. The noninteracting swimmer 
simulation shows the expected deviations from Cid at small r and large r due to the assumptions 
of the theory. For the interacting pushers, the tracer correlation is larger than that of the ideal 
theory and is system size-dependent as illustrated below, while for pullers, it remains close to that 
for independent swimmers, regardless of system size. The presence or absence of system size 
dependence in the velocity correlations is consistent with that for tracer diffusivity. 
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FIG. 3: (color online) The correlation function at = 3200 for different models at = 0.1: independent 
swimmers (blue), independent point dipole theory (black), full simulation of pushers (green), and pullers 
(orange), (inset) Correlation function for pushers at different system sizes with the spatial extent rescaled 
by k = i^-^'^L^-^^. The arrow denotes increasing system sizes = 400 — > 6400. 

The connection between extended spatial correlations and the system size dependent diffusivity 
for pushers is made clearer by returning to the system size dependent correlation length It. The 
inset to Fig. 3 shows the velocity correlation function for pushers at 0e = 0.1 for several different 
system sizes with distance scaled by It- Except at distances close to the total box size, this scaling 
collapses the results to a single curve. This suggests that the correlation length inferred from the 
diffusivity also controls the spatial velocity correlations. At intermediate distances, the curves are 
consistent with a generalized point dipole theory, with a correlation of the form C(r ) ~ / r 

The origin of the difference between pusher and puller dynamics is not completely clear, but a 
simple argument and some analysis can shed some light on the issue. Imagine an initially homo- 
geneous isotropic suspension of swimmers subject to a perturbation in the form of a shear flow. In 
a suspension of pullers, the shear flow distorts the orientation distribution and the resulting shear 
stress due to this distortion opposes the original shear flow, driving the system back to isotropy. 
This situation is closely analogous to a Brownian suspension of fibers or polymer molecules, the 
difference being that in the case of the swimmers the force dipoles that lead to the stress are in- 
trinsic rather than induced by the shear. In the case of pushers, however, the initial shear flow 
distorts the orientation distribution in the same way, but now the force dipoles have opposite signs 
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from the puller case and lead to a shear stress that enhances the original shear flow perturbation. 
This enhancement further increases the orientation of the swimmers thus implying instability with 
respect to shear perturbations of the homogeneous isotropic state. 

This argument can be made precise through a simple mean field theory (a related analysis 
without the above physical argument is given by Bill ). Let u = (v^) and ex. = (nn), where n is the 
director vector of a swimmer. For a homogeneous isotropic state, the number density c is constant, 
u = (n) = and a = In the point dipole limit the stress tensor generated by the swimmers 
is Tg = dcoL, where d is the dipole strength per swimmer (rf > for pullers, (i < for pushers). 
Stokes' equation relates u and Ts, and n evolves as would an infinitesimal material line subject to 
the constraint of constant length 132(1. Considering linear stability of the homogeneous isotropic 
state subject to a shear flow disturbance Uy{x) in the long wave limit. Stokes' equation and the 
evolution equation for axy become r]dxUy = —dca^y and dtotxy = {l/5)dxUy. Combining these 
shows that axy evolves exponentially in time with growth rate a = —dc/5r]; for pullers, the shear 
perturbation decays, while for pushers it grows, confirming the simple physical argument given 
above. This result is wavelength-independent in the long-wave limit and illustrates a mechanism 
for generating long-range correlations as seen in the full pusher simulations. 

The qualitative differences we observe between pushers and pullers naturally leads to the ques- 
tion of whether different organisms have evolved their method of swimming partially based on how 
significantly that method of swimming enhances transport in the fluid. In addition to biolo gica l 



35L 



371]. 



systems, a number of artificial micro-swimmers have recently been developed 1331 
The impact shown here of the mode of swimming on collective behavior could be an important 
design criterion for future devices of this kind. 
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